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Abstract. We address the issue of large-order expansions in strong-field QED. Our 
approach is based on the one-loop effective action encoded in the associated photon 
polarisation tensor. We concentrate on the simple case of crossed fields aiming at 
possible applications of high-power lasers to measure vacuum birefringence. A simple 
next-to-leading order derivative expansion reveals that the indices of refraction increase 
with frequency. This signals normal dispersion in the small-frequency regime where 
the derivative expansion makes sense. To gain information beyond that regime we 
determine the factorial growth of the derivative expansion coefficients evaluating the 
first 80 orders by means of computer algebra. From this we can infer a nonperturbative 
imaginary part for the indices of refraction indicating absorption (pair production) as 
soon as energy and intensity become (super)critical. These results compare favourably 
with an analytic evaluation of the polarisation tensor asymptotics. Kramers-Kronig 
relations finally allow for a nonperturbative definition of the real parts as well and 
show that absorption goes hand in hand with anomalous dispersion for sufficiently 
large frequencies and fields. 



PACS numbers: 12.20.-m, 42.50.Xa, 42.60.-v 



Large orders in strong-field QED 



2 



1. Introduction 

Recent years have seen a continuous progress in laser technology leading to ever 
increasing values of power and intensity. This is true for both optical lasers PjJ |2] 
and systems based on free electrons like the DESY vacuum ultraviolet free electron 
laser (VUV-FEL), a pilot system for the XFEL radiating in the X-ray regime [Sj. High 
intensities imply strong electromagnetic fields which approach magnitudes such that 
vacuum polarisation effects may no longer be ignored even at the comparatively low 
photon energies involved (1 eV ... 10 keV). The theory describing these effects is strong- 
field quantum electrodynamics (QED). 

The most exciting and therefore best studied effect arises when the ubiquitous 
virtual electron-positron pairs become real (Schwinger pair production [3]). This 
happens if the energy gain of an electron across a Compton wavelength X e equals its 
rest energy m e implying a critical electric field 

2 

777 

£ c =-^~1.3xl0 18 V/m. (1) 

This also is often named after Schwinger although it has first been obtained by Sauter 
upon solving the Dirac equation in a homogeneous electric field [3]. 

The physical situation to be analysed in this paper is as follows. We assume 
a background field consisting of an intense, focused laser beam of optical frequency 
fi ~ 1 eV <C m e . The associated gauge potential and field strengths are denoted 
and F^ u , respectively. The laser beam configuration is modeled by crossed fields with 
electric and magnetic fields constant, orthogonal and of the same magnitude, 

E = B = F = const , E JL B . (2) 

This configuration may be viewed as the zero- frequency limit of a plane wave, Q — > 0. 
In covariant notation one may write 

A^ = -^F^ v x , (3) 

with, for instance, F m = F = —F 3 i and all other entries vanishing. Working with 
crossed fields leads to enormous simplifications as will be seen below. 

The crossed-field configuration will be probed by a weak plane-wave field encoded 
in the potential a M and field strength = — d v a^. We denote the probe wave 
vector by 

k = u(l,nk), (4) 

where n > 1 represents the index of refraction (n = 1 in vacuo) and the unit vector 
k the direction of propagation. Note that n = 1/v with v = uj/\k\ < 1 being the 
phase velocity of the probe propagating through the modified vacuum. This set-up has 
recently been suggested for an experiment to measure vacuum birefringence for the first 
time [0]. 
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The corrections to pure Maxwell theory to leading order in the fluctuation a M are 
given by the effective action [7] 

SS = i J d 4 x d A y cv(x)n^(x, y; A) a v {y) . (5) 

This is basically determined by the polarisation tensor (or photon self-energy) IP" in 
the presence of the background field A^. 

In one-loop approximation the polarisation tensor is represented by the Feynman 
diagram 

11^ = wwQm 

where the heavy lines denote the fermion propagator dressed by the background field A 
(dashed lines below), 



S F [A] 



+ ... 



(7) 



Obviously, the first term on the right-hand side in (JJJ) corresponds to vanishing 
background fields, A = 0. In terms of the dressed propagator S F the polarisation 
tensor (in momentum space) is given by the integral 

U^(k;A) = -ie 2 tr J -0-^ S ' F (p)? 'S ' F (p - k) , (8) 

where 'tr' denotes the Dirac trace. 

The dressed propagator and the associated one-loop polarisation tensor are known 
exactly for a few special backgrounds only (see [H] and [§] for overviews). In particular, 
an exact solution is available for the crossed fields (J2J) due to Narozhnyi jTOj and 
Ritus ^T]. Their results may be written somewhat symbolically in terms of a spectral 
decomposition, 

W(k;A)= Mk; A) efer , (9) 

i=0,± 

with three orthogonal eigenvectors £j = ei(k; A) satisfying e» • k = 0. The eigenvalues 
Hi are given by somewhat complicated double integrals, the explicit form of which will 
be given below. It is important to note that the representation Q is valid for arbitrary 
external field strength and probe frequency. In other words, it contains all orders in the 
background intensity and probe energy. 

Let us briefly get an idea of the orders of magnitude involved. We assume that our 
probe is provided by another laser, possibly of high frequency, i.e. u ^> Q. For instance, 
if we think of the projected DESY XFEL as a probe with a projected wave length of 
A = 0.1 nm, hence u ~ 10 keV, we expect a ratio uo/m e ~ 0.02 1, that is low energy 
also for the probe. About the same frequencies can be achieved for photons from a 
laser-based Thomson back-scattering source |12j . 

Assuming the background as an optical Petawatt laser focused down to the 
diffraction limit the peak intensity will be EI] 

I = P/A 2 n ~ 3 x 10 22 W/cm 2 , (10) 
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for A = 2-k/VL ~ 1 /im and power P = 1 PW. Sauter's critical field (JTJ corresponds to 
an intensity of 

I c = E 2 C = m 4 Je 2 ~ 4 x 10 29 W/cm 2 , (11) 

still seven orders of magnitude larger than (fTU|) implying low intensity for our 
background laser J. Hence for present technology we have two primary small parameters 
[HI EI], namely 

v 2 = cu 2 /m 2 < 4 x 1(T 4 , (12) 
e 2 ee £ 2 /£ 2 =J/J c ~lx 1(T 7 . (13) 

We have listed the squared values as these typically arise as the leading-order (LO) 
contributions the reasons being essentially Lorentz and gauge invariance. 

In view of the small parameters ()12j) and (|T3*j) is seems safe to assume that LO 
accuracy in v 2 and e 2 will suffice for the time being. Anticipating further increase in 
laser power pQ we will, however, also consider higher orders in this paper and confront 
the outcome numerically with the exact (albeit implicit) one-loop results. This should 
provide some intuition about the limitations of derivative and weak-field expansions in 
cases where exact results are not available and one has to rely entirely on the accuracy 
of the expansions. 

The paper is organised follows. We first go through a general (mainly algebraic) 
discussion of the polarisation tensor where we also give the integral representations. 
Section El discusses the next-to-leading (NLO) results while Section |U extends these to 
order 80. Our findings are then compared with an asymptotic analysis of the integral 
representations in Sectional We finally present our conclusions in Sectional 



2. General Analysis of IP 1 ' 

In order to evaluate the following choice for the basis vectors 6j has been suggested 

in na, 

W EE F^ku , (14) 

If ee F^K , (15) 

d 1 ee F^by , (16) 

where, as usual, F^ u = (l/2)e tJ,ual3 F a p denotes the dual field strength. The first two 
basis vectors are obviously orthogonal to k, 

b-k = = b-k , (17) 

while the last one is not, 

c-k = -b 2 . (18) 

% While the presently attainable fields are 'weak' compared to the critical field they are certainly 
very strong by everyday standards. In particular they exceed static lab magnetic fields by orders of 
magnitude. This justifies the phrase 'strong-field QED' in the title. 
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(k fl k u \ 1 
9^ ~ -^r) c » = P ^ c - = -^ Pa kvb p b a , (19) 

the vectors b, b and d constitute an orthogonal dreibein which satisfies 
d^d v W b^b u 

-*+ir + ir = F ~- (20) 

The spectral decomposition Q may then be rewritten as 

jfiJv lulu im,im 

n- = n„^- + n + ^ + n_^, (21) 

and coincides with the one used in the monograph [8J. The eigenvalues IT depend 
on all the independent invariants one can form from the basis vectors and the 
background field strength. For crossed fields, however, most of these are either vanishing 
or dependent quantities. In particular, we have 

y = - = , (22) 

&=- \f^F^ = . (23) 

We are thus left with only two basic invariants, k 2 and b 2 , such that the eigenvalues in 
(|21|) can be written as 

n, = A; 2 P(fc 2 ,6 2 ) + ^(A; 2 ,6 2 ) , i = 0, ± , (24) 

where P and Pi are dimensionless polynomials in k 2 and b 2 (see below). Note that the 
&; 2 -term is the same for all eigenvalues. According to (jlj) we have 

k 2 = uo 2 (l-n 2 ) , (25) 

which is negative in a modified vacuum (n > 1). To calculate b 2 we note that the 
Maxwell energy-momentum tensor for the crossed background fields is 

= F^ X F X U - g^y = F imX F x u , (26) 

whereupon b 2 simply becomes 

b 2 (k) = -T^k^K . (27) 
Using (@J) again this translates into the expression 

b 2 = -uj 2 [H(l + n 2 ) -2nS-k-n 2 H k ] , (28) 
where we have introduced the abbreviations 

H =T 00 = \{E 2 + B 2 ) = F 2 , (29) 

S i _ T oi = e m E j B k ^ ^ 

H k = (E. k) 2 + (B ■ k) 2 , (31) 
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the first two of which represent the Maxwell energy- density and the Poynting vector, 
respectively. Note that both these quantities are directly related to the intensity, 

H=\S\ = I. (32) 

Expression ()28|) for b 2 coincides with the quantity on p. 21 of |Sj. There is, however, an 
alternative expression which nicely illustrates the kinematics involved. Introducing the 
background 4- vector K = 0(1, K), K 2 = 1 and assuming background gauge d ■ A = 
one finds 

b 2 = b 2 (6) = -oo 2 I(l -ncosO) 2 < , (33) 

where 9 is the angle between the probe and background directions (k and K, 
respectively) . Note that b 2 becomes independent of n for a perpendicular configuration 
implying b 2 (%/2) = —uj 2 I. 

In order to obtain reasonably simple expressions for the coefficient functions P and 
Pi multiplying k 2 and b 2 /I c in ()24)1 we follow ^JJ and trade the invariants k 2 and b 2 for 
the dimensionless parameters 

A = k 2 /m 2 e = v 2 (1 - n 2 ) < , (34) 
k 2 = - b 2 /I c m 2 e = e V (1 - n cos 6) 2 > . (35) 

The eigenvalues of IP" from (|2~4~j) are then given by the expressions 

n (A,K 2 ) = m 2 e XP(X,K 2 ) , (36) 
n±(A, k 2 ) = n (A, k 2 ) - m\ k 2 P±(A, k 2 ) , (37) 

implying P = in particular. Using the integral representations given in [IT] the 
remaining coefficient functions P and P± may be compactly written as 

2a f 1 

P{X,k 2 ) =P(A,0) / dxxxf^z), (38) 



7T 



2a f 1 

P(A,0) = — / dxxx ln(l- \xx) , (39) 



7T 



/'W • (40) 

In the above, we have introduced the Feynman parameters§ x and x = 1 — x and the 
variable [TT] 

z = z(A,k,x) = 1 zrfS , (41) 
(kxx) 2 /' 3 

which is the argument of the auxiliary functions 

POO 

f(z) = % \ dte- izt ~ it3/3 = 7i [Gi(z) + iM{z)\ , (42) 



o 



>3C 



/,.(:) ^ / ^ e"** ( C -*V» - 1 ) . (13) 



o 



§ The Feynman parameter x is related to the integration variables v in jH] and i> in via v = 2x — 1 
and i; = 1/xx, respectively. In particular, dv = (x — x)/(xx 2 ) dx. 
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/' denotes the derivative of / with respect to z; Gi and Ai are Scorer and Airy functions, 
respectively (see e.g. JE] or the Digital Library of Mathematical Functions [T7j). 

As both Gi and Ai are real for real arguments we can immediately infer from (|4()jl 
and (|4"2j) that the eigenvalues will develop imaginary parts determined by the Airy 
function Ai(z) in the integrand. It is thus natural to expect that also the index n of 
refraction will become complex, its imaginary part signalling absorption of photons by 
the vacuum. Of course, the only physical interpretation of this phenomenon is pair 
production. 

Enormous simplifications arise naturally for vanishing external field, i.e. k 2 = 0. 
In this case all eigenvalues become degenerate, IT = II = m 2 . XP and one just obtains 
Schwinger's one-loop expression for the photon self-energy [3] in the form 

W u = U{k 2 ) P"" = m\X P(A, 0) P^ . (44) 

The scalar function P(X, 0) is given in (ffilj) and coincides with standard text book results 
obtained via covariant perturbation theory (see e.g. jTH], Ch. 7). For small momenta it 
becomes 

P(A,0) = -— A + 0(A 2 ), (45) 

107T 

so that I1(A; 2 ) = 0(k 4 ) for vanishing external field. 

The next logical step is to expand the exponentials in (|4~2jl and (}4*3|) in powers 
of A and k 2 . This amounts to a derivative expansion in the probe field (keeping the 
background fixed) and will be further pursued in the next section. 

Before that let us conclude the general reasoning by determining the dispersion 
relations for k. Adopting a plane wave ansatz for the probe field, f^ u = i{k^e u — 
k^e^) exp(zfc • x), and adding the Maxwell term yields the wave equation 

nTe v = -{k 2 ?^ - n^)e y = , (46) 
which has nontrivial solutions only if 

detD^(fc) = 0, (47) 
i.e. if an eigenvalue of vanishes, 

k 2 - IL,(A; 2 , b 2 ) = h^(k 2 } b 2 )k^k v = . (48) 

Following ^1 we have introduce effective metrics to make explicit that (|4*Hj) 
represents 'modified light-cone conditions' jHj. Inserting (|24|) the three metrics become 

hT(X, k 2 ) = <T [1 + P(A, k 2 )] + T^Pi(X, k 2 )/I c . (49) 

As Pq = the metric is conformally flat, 

h^(X,K 2 ) = g^[l + P(X,K 2 )} , (50) 

and does not modify the light-cone. Note in particular that in contrast to a plasma 
of real charges weak crossed fields do not generate a longitudinal photon so that 
remains unphysical. 
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The other two metrics, however, are physical and nontrivial leading to a modified 
light propagation. As h\ 7^ h 2 this implies 'birefringence of the vacuum'. For vanishing 
external fields (k 2 = 0) all metrics merge into the conformal metric Hq U {\, 0) which 
describes a standard light-cone. 



3. NLO Derivative Expansion for the Probe 



The integral representations (JH8H40J) for the polarisation tensor may be expanded in 
powers of A and k 2 which are both 0(v 2 ). This derivative expansion becomes particularly 
straightforward if we first expand in v and only afterwards perform the integrals, a 
procedure adopted throughout Sections 01 and |H Thus, we first rewrite the derivative of 
and the second function (}4*3*j) in order to exhibit the dependence on k and A, 



f'(z) = (kxx) 4/3 / dre iXxST e- tT - iK2xW/3 , (51) 



h{z)= J ^e iA ^ T e- iT ^e- iKWr3/3 -l) . (52) 

Upon expanding the exponentials the integrals both over Feynman parameters and 
proper time become elementary so that the eigenvalues Ilj to 0(z/ 4 ) are found to be 

U /ml = - t^A 2 + -£-\k 2 , (53) 

157T 1057T 

U + /ml = -^. 2 - -^-A 2 - ^X. 2 - , (54) 

+/ e 45vr 15vr 315vr 945tt ' v ; 

TT / 2 ^tt 2 a \2 ^ a \ 2 ^ a 4 /rr\ 

Tl-/mz = k A \k k . 55 

1 e 457T 157T 45vr 135tt v ; 

Note that only Il± have LO contributions, Il± = 0(n 2 ) = 0(u 2 ), whereas n = 0(A 2 ) = 

0(z/ 4 ). To proceed we specialise to a head-on collision for which birefringence becomes 

maximal, cf. (}3*3*|) and (|33|) . In this case, k ■ K — — 1 so that k 2 attains its maximum 

value, 

K 2 = e V(l + n) 2 . (56) 

The nontrivial dispersion relations \3± = A: 2 — n=b = implicitly determine two nontrivial 
indices of refraction n± as functions of the small parameters e 2 and v 2 . To solve the 
equations D±(n, e, v) — for the indices n± we write 

n± = 1 + A± = 1 + e 2 5± . (57) 

457T 

The last expression takes into account that corrections to n = 1 are due to nonvanishing 
external field intensity e 2 the coupling to which is 0(a). The remainders 5± are expected 
to be of order unity. For the value e 2 = 10 -7 from (|13|) the prefactor is of order 10~ n 
so that at present the deviations from n = 1 are extremely small. It hence remains an 
experimental challenge to really measure them |d1 I20j. 
If we expand the deviations according to 

6± = 5 0± + 5 2± z/ 2 + O(z/ 4 ) . (58) 
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we find at LO in k 2 or v 2 i.e. 0{u° 

14 . f <\( 



2 



14<M +7^ + 0(a 2 e 4 ) \ , (59) 



+ l-7ae 2 /457r [ 45vr 

^- l-4lv45^ 8 { 1+4 ^ + Q(a2e4) }- <60) 
Note that the first terms on the right-hand side are exact to all orders in the intensity 
e 2 . An independent check of the LO results based on the Heisenberg-Euler Lagrangian 
will be performed in App. A. 

To the best of our knowledge, the coefficients of order e° in (foT?|) and (f6TJ|) have first 
been obtained in Toll's thesis [21] || and independently in [23 121] (see also [7] ITU] 123]). 
The NLO=0(z/ 2 ) expressions are somewhat more complicated, 

'416 184a 392a 2 \ e 2 

(61) 



21 45vr 675tt 2 J (1 - 7ae 2 /45vr) 4 ' 
32 64a 128a 2 \ e 2 

T ~ 45tt + 675tt 2 J (1 - 4ae 2 /45vr) 4 ' ^ ^ 

Again, these expressions are exact to all orders in e 2 . For later purposes it is useful to 
attach names to the different factors, 

h± = c 2 ±{a)e 2 s 2 ±{ae 2 ) . (63) 

According to ()61|) and (j62j) the c 2 ± are polynomials in a/457r while the s 2 ± have the 
series expansions 

28a 2 98a 2 4 fi . . 
s 2+ = 1 + — e 2 + — ^e 4 + O e 6 , 64 

16a 2 32a 2 4 _ . fi . . 
s 2 _ = 1 + — e 2 + — ^e 4 + O e 6 . 65 
457r 405^^ 

Note that both at LO and NLO the QED expansion parameter is actually a/457r ~ 
5 x 10~ 5 rather than a itself. The coefficients c 2 ± are dominated by the leading terms 
of 0(a°) which both have a positive sign. This implies that (at least to first nontrivial 
order in the derivative expansion) the indices of refraction increase with frequency so 
that we have normal (rather than anomalous) dispersion. 

Summarising the findings above we see the following pattern emerging. Each of the 
two indices of refraction has a derivative expansion in z/ 2 , 

2 00 

n± = 1 + A± = 1 + ^- V 5 2l± (a, e) v 21 , (66) 

457T L — ' 
1=0 

where the S 2 i can be factorised by generalising (J63|) . 

5 2l± = c 2l ±{a) e 2l s 2l± {ae 2 /A5n) . (67) 

|| As we have not been able to get hold of this unpublished work our statement is based on the notes 
[22] where the relevant figure of Toll's thesis is reproduced on p. 33. 
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Hence, the different coefficient functions have the generic behaviour 

S21 = 0(e 21 ) , (68) 
c 2l = 1 + 0(a) , (69) 
s 2l = 1 + 0(ae 2 ) , (70) 

for both subscripts ±. Keeping only the leading orders in a in (157)1 . i.e. 5 2 i — c 2 i(0)e 21 , 
results in a compact expression for the NLO derivative expansion, 

5± = 11 ± 3 + 32 ° ^ 96 eV + 0(a) + 0(ae 2 ) . (71) 

Again, the relative plus sign between the first and second terms signals normal 
dispersion. Extrapolating (|7T)l to all orders adopting the same approximations we expect 
(158]) to become a function of ev only, 

00 

<^ J>(0)(e*,f. (72) 
2=0 

If we have a closer look at the coefficients co± and c 2 ± (at 0(a°)) we see that they 
actually increase, the ratios being almost the same for both indices, 
c 2+ 416/21 c 2 _ 32/3 

~ ~ 1.4 and ~ ~ 1.3 . (73) 

c 0+ 14 c - 8 

This is a first hint that our expansion in v (or ev) is only asymptotic which has to be 
expected upon comparing with closely related derivative expansions of effective actions 
O EE] • In the remainder of the paper we will investigate this issue in detail. 



4. Large-Order Derivative Expansion for the Probe 

Given the present day power of computer algebra systems we have decided to actually 
check whether normal dispersion persists beyond NLO in the derivative expansion. The 
answer is affirmative as shown in the Mathematica plots of figures [I] and 121 They display 
the second, sixth and tenth order in the derivative expansion of A = n — 1 as a function 
of v at fixed background intensity e 2 = 0.1. The LO results ()59j) and flOOl) are given by 
the common intercept of the different graphs. One clearly identifies normal dispersion 
and a tendency for the curves to diverge with increasing order. 

As already stated the results of the preceding section indicate strongly that the 
expansion of the index of refraction is asymptotic both in frequency and intensity. 
Hence, for both the derivative and weak (external) field expansions we expect factorial 
growth of the expansion coefficients at large orders. 

To test this expectation numerically we expand A± in accordance with (jfifij) . 

A ± (^,e) = ^J> ± (^, (74) 

l>0 
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0.2 0.4 0.6 0.8 1 V 

Figure 1. Higher order results for the deviation A + of the index n+ from unity as a 
function of v — aj/m e at hxed background intensity e 2 = 0.1. Full line: second order; 
short-dashed line: sixth order; long-dashed line: tenth order. 




0.2 0.4 0.6 0. 
Figure 2. Same as Fig.nfor A_. 



1 V 



with S 2 i± as given in (J67j) and form the ratio^f 

$2l+2,± C 2 Z+2,±(0) r 
T2l± = — = J-T- [1 



+ 0(a) + 0(at 



(75) 



e 2 S 2i ± c 2 ;,±(0) 

The coefficients can (in principle) be determined by expanding the integral 
representations (jHBHUll) using the auxiliary functions (J4"2j) and (f4"3j) (see next section). 
For simple factorial growth, 621 ~ r(/), r^i would depend linearly on I. 

In what follows we want to ask the question what can be said about the factorial 
growth without using any knowledge of the special functions (|4"2j) and (|4"Hj) . That is, 
we perform a brute-force derivative expansion before calculating any integrals, trying 
to extend the method of the previous section to orders as large as possible to gain a 
maximum of information. The philosophy behind is our expectation that for realistic 
backgrounds such as laser fields this type of derivative expansion will be one of the 
main tools when analytical methods are not available. We thus expect our method 

We thank Paul Rakow for bringing this idea to our attention [37] . 
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to provide results complementary to numerical approaches like world-line Monte Carlo 

[2E1 123 EDI Ell- 
in order to extract the asymptotic behaviour of 621 or t2i with high accuracy it is 
obviously desirable to go to the largest orders possible. Again, this is only feasible with 
the aid of computer algebra and was performed with an optimised Mathematica routine. 
Using a standard desktop PC we were able to achieve a maximal order of i = 82. The 
determination of all 82/2 coefficients for a given value e 2 of intensity takes about an 
hour. 

In Fig. s El and |U we have plotted the ratio ([75]) as a function of order 21 for different 
values of the intensity e 2 ranging from 0.1 (subcritical) to 500 (supercritical). 




Figure 3. Successive coefficient ratio r2i+ as a function of order 21 for four different 
intensities, e 2 = 0.1 (♦), e 2 = 1 (*), e 2 = 100 (■) and e 2 = 500 (A). 




Figure 4. Successive coefficient ratio r^i- as a function of order 21 for four different 
intensities, e 2 = 0.1 (♦), e 2 = 1 (*), e 2 = 100 (■) and e 2 = 500 (A). 

One notes the following: First, the e-dependence of r 2 i is rather weak as expected 
from the discussion in the preceding section and the second expression in (|75j) . To 
recognize corrections to the leading behavior which is independent of e one has to go to 
very large intensities and/or orders. For instance, the graphs for e 2 = 0.1 and e 2 = 1 
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basically are on top of each other. One needs e = O(10 2 ) to see 10% deviations from the 
graphs with small e values. This is consistent with the coefficients of e 2 in (|75j) being of 
order 10" 3 , cf. fgj) and ffgfl . 

Second, and more importantly, the ratio r 2 ; seems to depend quadratically on I 
which is simply explained by an asymptotic behaviour r 2 \ ~ T(2l). Let us try to perform 
a more quantitative analysis. The following general ansatz 



c 2l (Q) = p M Y(2l - a) + tT(21 - a - I) 



(76) 



which is tailored after examples where exact results are available {26] E2] , implies a ratio 



r 2 i = P 



(2l) 2 + (l-2o) 21 + a 2 — a — 2t\ = a>2 (2l) 2 + ai 21 + a . (77) 



Fitting the curves of Fig.sElandE]to quadratic polynomials we find the a$ values listed 
in Tabled 



Table 1. Fitted values for the polynomial coefficients in J77J). 





e 2 = 0.1 


e 2 = l 


e 2 = 100 


e 2 = 500 




0.56256(11) 


0.56295(35) 


0.607(11) 


0.8383(38) 


a 1+ 


1.6892(29) 


1.6899(49) 


1.817(44) 


2.51(42) 


a 0+ 


1.691(41) 


1.693(46) 


1.69(27) 


1.4(2.5) 


a 2 - 


0.562535(27) 


0.56275(15) 


0.5865(15) 


0.699(11) 


ai_ 


1.68741(46) 


1.68795(16) 


1.7591(21) 


2.10(15) 


a _ 


1.3926(36) 


1.3928(70) 


1.38(13) 


1.25(90) 



The errors have been estimated by also fitting higher-order polynomials and 
monitoring the change in the coefficients. For all fits an extrapolation to / — > oo has 
been performed using Laurent series techniques. 

Via (|77jl the values of Table [T] translate into the asymptotic coefficients p, o and r 
introduced in (fT6|) as listed in Table |2l 

Table 2. Numerical values for the polynomial coefficients in J2jJ. 





e 2 = 0.1 


e 2 = l 


e 2 = 100 


e 2 = 500 


p+ 


0.75004(15) 


0.75030(47) 


0.779(14) 


0.9156(42) 




-1.0014(23) 


-1.009(34) 


-0.9967(91) 


-1.00(24) 




-0.501(65) 


-0.502(69) 


-0.40(35) 


0.2(2.2) 


P- 


0.750023(36) 


0.75017(20) 


0.7658(20) 


0.836(13) 


(7_ 


-0.99983(34) 


-0.99973(26) 


-0.9997(20) 


-1.002(84) 


T_ 


-0.2380(52) 


-0.238(12) 


-0.18(22) 


0.11(99) 
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One notes that p and a are the same for both indices of reflection i.e. p + = p_ = p 
and a + = o~_ = o~ while r + 7^ r_. Furthermore, <r seems to be independent of intensity 
(unlike p and r). We can even guess the following 'analytic' values for p and a, 

p = 3/4 , and o = -1 . (78) 

For e 2 < 1 these values for p and cr should be quite accurate. Not unexpectedly, the 
largest errors reside in the sub leading coefficients t±. In view of this we refrain from 
any further guesses and just quote the r- values of Table El 

With the coefficients being determined let us plug the ansatz (|7U|) into (|T4*j) adopting 
the approximation (jfi7j) so that the expansion parameter becomes ez/ rather than v. This 
should be fine for e 2 < 1. As our ansatz (|7B|) is supposed to hold only asymptotically 
for large orders we cannot expect to describe the low orders accurately. We nevertheless 
proceed by eliminating the Gamma functions using the integral expression 

r°° e -^/ps 

p z T{z)= / ds^^. (79) 
Jo s 

This transforms the required summations over I into geometric series and yields the 
following integral representation for (|74|) . 

2 /*oo 

A±(e, ") = N ± %L ds - S [(psT + T ± (psy +1 ] e- 1 ^ , (80) 

where we have allowed for an undetermined scale N±. As it stands the integral is 
ambiguous as there are poles on the real axis. The left-hand side may be viewed 
as originating from the derivative (or tv) expansion which has real coefficients. It is 
therefore tempting to interpret (|8()|) as the following dispersion integral (see e.g. [33J, 
Ch. 18), 

ReA ± (ez/) = N±-% / ds - . . 9 ImA ± (s) , (81) 



7T ' Jo s 2 — (ez/) 2 

where the pole ambiguity has been resolved in terms of the Cauchy principal value 

denoted by ^2. The right-hand side contains a 'nonperturbative imaginary part' [3*2*] . 

2 

Im A ± (ez/) = N ± ^- {evpf (1 + r ± ez/p) e^ 1 ^ , (82) 

i.e. an exponential that cannot be obtained from a perturbative expansion in ez/. 
Inserting the universal values (|78J) and r± from Table El this becomes 

a / \ , T ae 2 4 / 3ez/ f 0.501 \ 4 , w , , 

lmA ^ = ^90 to I 1 " X {0.24}) 6 ' ' < 83 > 

Note that the small parameter involved seems to be 3ez//4. 

Summarising we can say that our large-order derivative expansion has provided us 
with a quantitative formula for the asymptotics of the series which by means of (JHTj) 
could be turned into a statement about the nonperturbative imaginary part. As already 
stated, this should be directly related to absorption and pair production. In what follows 
we will check our findings by an analytic discussion of the integral representations of 
the eigenvalues IT, in particular of (jJOJ). 
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5. Analytic Results 

Proceeding analytically is equivalent to analysing the (derivative of the) auxiliary 
function f(z) introduced in (}4*2*|) . The leading orders for both weak and strong external 
fields have already been investigated by Narozhnyi (10). In this section we want to go 
further and discuss arbitrary large orders in the low-frequency, weak-field expansion 
of the eigenvalues ITj. Nevertheless, to keep things manageable we follow Narozhnyi's 
example and adopt an approximation that leads to substantial simplifications. Namely, 
we neglect vacuum modifications in expressions involving the eigenvalues Hi(\, k) by 
setting n = 1 or, equivalently, A = in their arguments. From (jHBj) and (|37j) we thus 
have the approximate identities 

it (a,k 2 ) ~n (o,/€ 2 ) = o, (84) 

n ± (A, k 2 ) ~ n±(0, n 2 ) = -m 2 e K 2 P±(0, k 2 ) . (85) 

According to (jHfij) we have in the same approximation, 

n — lev , (86) 

which turns our eigenvalues into functions of the product ev in line with our discussion of 
the previous section. The upshot of all this is that the dispersion relations ()48|) simplify 
drastically so that one is left with the task to determine P± in 

A±(k, A) ~ A±(ei/) ~ 2e 2 P ± (ev) . (87) 

As an aside we remark that according to (J3*3j) . (|34|) and (|35|) Narozhnyi's approximation 
is expected to be particularly good if probe and background are perpendicular as b 2 and 
hence k are indeed independent of n in this case" 1 ". The perpendicular setup has recently 
been suggested as a means to look for 'vacuum diffraction' [31] . 

To evaluate (|57jl we want to calculate P± as a power series in ev using the integral 
representation (|40jh From the definition of k 2 in (}35|) or ()86j) and of z(A, k,xx) in (j4T| 
it is clear that an expansion in ev corresponds to determining the large- z asymptotics 
of 

f'{z) = 7rGi 7 (z) +iitki'{z) , (88) 

which appears in the integrand of (|40|). Using asymptotic expansions for Ai(z) given in 
[in] and Gi(z) in jHSl one obtains for (fKKj) 

1 1 oo . oo 

/'w = 4 - 7 E ^ - \ D- 1 )*^ c-* , (89) 

fc=0 fc=0 

with the abbreviations 

C =^ 3/2 , (90) 

= (3^ + 4^(3^ + 3) (91) 

" 3 fc r(fc + i) ' 1 } 

_ Qk + i r(3fc + i/2) 

fc 6ife-l 54 fc fc! T(k + 1/2) ' 1 1 

+ with the replacement 2eu — > understood in (|86|) 
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The expression flH9~|) now has to be plugged in the integral (J4TI|) . Throughout the 
subsequent calculations we use Narozhnyi's approximation which implies 

z ~ (kxx)- 2/3 , C ^ 2/3kxx . (93) 

In what follows we will separately calculate real and imaginary parts of P±. 

5.1. Calculation of the Real Part 

Using (JHEJ), ((SS| and the integral representation of the Beta function, 

1 dx x™- 1 x m - 1 = B(n, m) = T ^ T{ - m ] (94) 

the real part of P± can be obtained in closed form, 

2 

ReP±(/c) = ^^G^ 2/ {2E(2/ + 2,2/ + 2) + (l±3)5(2/ + 3,2/ + 3)} . (95) 
z>o 

Here, we have introduced the new expansion coefficients 

a "{ 1, 'ill (96) 

^ gi-i , I > 

The leading term (/ = 0) in the expansion (J9*5j) yields for (JKTj) 

A ± (ez/)=2e 2 P ± M^(ll±3)— , (97) 

which reassuringly is consistent with (|71|). 

It is the presence of the Beta functions in (J9*5j) that causes factorial growth. To 
analyse this in the spirit of the previous section we expand the corrections to the indices 
of refraction according to (JB5J) and (J72J) . 

Re A± = 2e 2 ReP ± = ^ ^ c 2Z , ± (0) (e^) 2i . (98) 

«>o 

Comparing with (|95p we find the coefficient ratios 

_ 2P(2/ + 4, 2/ + 4) + (1 ± 3)B(2l + 5,21 + 5) 

r2,± ~ G, 25(2/ + 2,2/ + 2) + (1 ±3)5(2/ + 3,2/ + 3) ' ^ ' 

These can be evaluated exactly to give 

2(6/ + 13)(3/ + 4) (31 + 2) (21 + 3)(/ + 1) 



(6/ + 7)(4/ + 9)(4/ + 7) 

= A(2/) 2 + ^2/ + ^-^I + 0(l// 2 ), (100) 
16 V ; 32 256 256 2/ y 1 ' ' K ' 

(6/ + 14)(3/ + 2)(2/ + 3)(/ + l) 
r-2i- = 



(4/ + 9)(4/ + 7) 

A ( 2 /)2 + ^1 2/ + — - — I + 0(1// 2 ) . (101) 
16 V ; 32 256 256 2/ v 1 ' K ' 
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Obviously, as the r 2 i± are independent of intensity they cannot describe all curves of Fig.s 
El and 0] which, after all, are intensity dependent. However, we expect good agreement 
between numerical and analytical values of r 2 i± for small intensities e 2 < 1. This is 
nicely corroborated by the graphs of Fig. El where the results for e 2 = 0.1 are on top of 
each other for all I. This near-perfect agreement is due to the fact that upon neglecting 
the e dependence of the coefficients (16 7|) (and within Narozhnyi's approximation) all 
expansion coefficients are basically factorials, cf. (jHU). 



5000- 



4000- 




10 20 30 40 



/ 

Figure 5. Coefficient ratio T2i+ as a function of order I. Diamonds: numerical values 
for e 2 = 0.1; crosses: numerical values for e 2 = 500 (both as in Fig. |3J); full line: 
analytical ratio from l|99() . 



Matching (jlOOJ) and (jlOlj) with (J77)l we read off the coefficients 

a 2 = — = 0.5625 , (102) 
16 

ai =^ = 1.6875, (103) 

452 , . 

a 0+ = — = 1.765625 , 104 
256 

356 . , 

a o- = ^ = 1-390625 , 105 
256 

which compare favourably with the entries of Table^for e 2 < 1. The cij translate exactly 

into the asymptotic parameters p = 3/4 and a = — 1 we had guessed in (|7S|) while the 
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r-parameters become 

r + = - !| = -0.5694 , (106) 

r = - H = -0.2361 . (107) 

Again, these agree satisfactorily with the entries of Table |2]for small e, in particular for 
subscript 'minus'. 

It remains to rewrite our nonperturbative imaginary part (J83|) by replacing the 
numerical r± by their analytic counterparts just obtained, 

_ a / \ 1T «e 2 4 / 1 3ez/29±l2\ 4/3ei , , ino , 

ImA±(«/) =N± l e~ A '^ u . (108) 

V ; 90 2>ev V 4 72 J K ' 

Let us finally try to check this by a direct calculation of the imaginary part. 
5.2. Calculation of the Imaginary Part 

The determination of the imaginary part of P± is more involved as the integrand ()40|) 
contains an exponential that cannot be expanded in powers of xx. Explicitly, the 
integrand is of the form 

e~ c {xx) n ~ exp(-2/3rax) {xx) n , (109) 

which is to be integrated from to 1. In order to proceed analytically we note that the 
product xx is peaked at the value 1/4. Thus, it seems feasible to perform the integral 
via saddle point approximation which results in the formula 

1 / _ \ 1/2 

-axx{^, — \b / \ A—b^—Aa 



I(a,b) = / dxe- axx (xxf ~ 4" b e- 4a . (110) 

Jo \l<oa + AbJ 

Here, the first prefactor has been obtained by extending the integration to the whole 
real axis in order to have a Gaussian integral. In our case, the parameters a and b are 
given by 

a = 2/3k ~ l/3ei/ and b = k ± 1/2 , k = 0, 1, 2, . . . . (Ill) 

Comparison with a numerical evaluation for k = ... 5 shows that the error of the 
approximation is about 10% for a = 1, 1% for a = 5 and 0.1% for a = 10. Hence, for 
large a (which we have) the formula (jllOj) should work very well. 

Having gained sufficient confidence in our saddle point approximation we rewrite 
the imaginary part of (jjOJ) using ([89)1 . 

?t ~ ^ *T 3 / 2 ^(-l) fc A {2 /(a, fc - 1/2) + (1 ± 3) /(a, + 1/2)} , (112) 

fc>0 

with Ak as defined in (|92j) . Employing our integration formula (jllOj) we find 

^M&^Z^(- 3 i)*' (113) 

fc>0 v 7 
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where the new expansion coefficients are given by the somewhat lengthy expression 

(114) 



„ _ . . 9 ±3 72k 



9 ± 3 -7±3 
4 ' (6Jfe-7)(6ife + l) ~^ {> '~ ' ' ; 

Admittedly, this is not too illuminating so we evaluate the series (J113)) to order k 2 or 
(ez/) 2 and use (|H7|) to determine 

T a ir «e 2 4 4/ ~ f 25^123ez/ 265 =f 168 /3ez/ \ 2 1 , iir . 

ImA ± ~iV ±sp e _4/3eI/ <l + — - + ^zr-z-l ■ ■}. (115) 

' p 90 3ez/ [ 72 4 2 x 72 2 \ 4 / J v 7 

In the above we have introduced the saddle point normalisation factor* 

tf ± „ = (3±l)gy|. (116) 

We refrain from identifying this with the analogous factor in ()108j) as the approximation 
schemes involved do not allow for a direct comparison of the normalisation. The latter 
will be fixed in a moment by means of the dispersion relation (|8ip. 

Putting these niceties aside we stress that we find perfect agreement between ()108j) 
and ()115j) concerning the LO dependence on powers of eu and the nonperturbative 
exponential. The only discrepancy resides in the r-coefficients multiplying the sub- 
leading terms of order eu. In view of the fairly different approximations employed this 
is probably not too surprising. 

6. Discussion and Conclusion 

With the consequences of our large-order analysis confirmed analytically we can go even 
one step further. We can use the Kramers- Kronig relation (|8*T|) to actually define the real 
part of the indices of refraction nonperturbatively, i.e. without relying on the derivative 
expansion. To this end we take the leading order of the imaginary part ()108j) or ()115j) 
as a model by writing 

2 A 

Im A±(e, u) = N ± — e" 4 ^ . (117) 

Plugging this into the the dispersion relation (}8T| we can actually perform the principal 
value integral analytically with the result 

90vr 3eu I 

Matching the perturbative small-ez/ behaviour fixes the normalisation to be 

iV± = ll±3. (119) 

Interestingly, in (|118J) exponential integrals Ei appear multiplied by exponentials, 
constituting a paradigm example of functions displaying factorial growth expansion 

* The leading term in l|115|) has already been determined by Ritus |1 lj . His normalisation differs from 
ours by a factor of 2, a discrepancy we have not been able to trace. The most plausible explanation 
seem to be the ambiguities in determining the prefactor in the saddle point approximation. 



Re A± = N± —— — |Ei(4/3ei/) e~ 4/3 ^ - Ei(-4/3ei/) e 4/3ei 4 
907T 3ev I J 



:ii8i 
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8e±05 - 



6e±05 

A + 

4e±05 



2e±05 - 




Figure 6. Real and imaginary part of the deviation A + of the index n + from unity 
(dashed and lower full curve, respectively). The upper full curve bending upwards is 
the series expansion of A + to second order (same curve as in Fig^). 



coefficients. With the real part thus determined we summarise our findings in Fig. |U] 
which shows both real and imaginary part of A + as a function of frequency v for fixed 
e 2 = 0.1. We have added the series expansion to second order (the full line of FigJTJ) 
which coincides well with the exact real part for small v where the factorial growth is 
not visible yet. 

The graph for A_ looks almost identical with a slight shift in the vertical scale due 
to the difference in normalisation. Hence, we refrain from producing an extra plot. 

By construction, the real and imaginary parts of A (thus also of n = 1 + A) are 
related by the Kramers- Kronig relation (|81jl with Imn ^ signalling absorption, i.e. 
pair production. By looking at Fig. |H]we see that the latter sets in roughly at the critical 
value of v = 1. Here, Ren attains its maximum and decreases for v > 1 while Imn 
increases. Mathematically, we may state 



<9Ren-i- 

> for uj <C m e , (normal dispersion) 



duj 
<9Ren ± 
dui 



(120) 

< for uj > m e , (anomalous dispersion) 
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where we have reinstated physical units (recall that u is the probe frequency). We 
conclude that it is a consequence of the Kramers-Kronig relations based on the 
fundamental principle of causality that absorption is intimately connected to anomalous 
dispersion. 

The physics involved can be wrapped up as follows. We have analysed the influence 
of crossed background fields on the propagation of light (e.g. laser beams). Using exact 
integral representations for the eigenvalues of the polarisation tensor we found that to 
all orders in probe frequency and background intensity there is birefringence of the 
vacuum induced by the crossed background fields. The effect can be described in terms 
of background dependent effective metrics h± implying dispersion relations h^k^k u = 
which describe distorted light-cones. Solving for the indices of refraction n± one finds 
that they are frequency dependent, starting out with normal dispersion. At critical 
energy and intensity anomalous dispersion sets in together with absorption due to pair 
production. 

At present an experiment is being designed that plans to measure vacuum 
birefringence using a high power laser background probed by x-ray beams j^j. The 
experiment is quite demanding as one has to measure ellipticity signals with a sensitivity 
at the order of 10~ n for presently available probes and lasers. As the signal is 
proportional to u 2 e A it may be readily enhanced by increasing probe frequency and 
background intensity. Within the next few years one expects a reduction of the required 
sensitivity down to 10~ 4 . Normal dispersion is an NLO effect, hence implies a signal 
of order z/ 4 e 4 requiring a sensitivity of 10~ 8 within the envisaged scenario. This is still 
within the theoretical limits of measurability [36J. Pair production shows up in terms of 
a nonvanishing imaginary part or as anomalous dispersion of the real part. To become 
observable both effects require parameters v and e close to their critical values. For 
instance, if e ~ 10~ 2 then ImA ~ 10~ n for v ~ 0.8. Even these moderate values cannot 
be attained at present so that it is presumably more reasonable to look for positrons 
rather than optical signals to detect pair production. 

From a theorist's point of view one should also consider two-loop corrections and 
the influence of nonconstant backgrounds. The two-loop corrections to the Heisenberg- 
Euler Lagrangian (cf. App. A) have been calculated by Ritus jHZj- They amount to a 
replacement of the LO coefficients of (J7T|) according to 

f 14 If 14 (1 + 13150/252^)1 

\ 8 J \ 8(1+ 40a/9vr) J ' 1 ' 

These are both one-percent corrections to the LO=0(v°) behaviour in (|71jh 

Regarding nonconstant backgrounds it is possible to slightly relax the crossed-field 
assumption in a perfectly controlled manner. Laser beams may be more realistically 
described as Gaussian beams rather than plane waves. The former have a Gaussian 
profile in transverse direction of 'waist size' Wq and a Lorenz profile in the longitudinal 
direction characterized by the 'Rayleigh length' zq. One can form the small dimensionless 
parameter A = wq/2z <C 1 [38J which describes the deviation from the crossed-field 
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limit corresponding to A = 0. Naturally, one expects that there will be 0(A)-corrections 
to the results presented here. 

It is this context of nonconstant backgrounds where the intuition gained in this 
paper is expected to pay off. In this more realistic case, we cannot hope to have any 
exact analytical results available. Thus it is important to know both the region of 
validity and the limitations of derivative and weak-field expansions. From our results 
they are both expected to break down if the product of frequency squared and intensity 
becomes u 2 I = 0(m®/e 2 ), or, in dimensionless units, e 2 u 2 = 0(1) (see Fig. EJ. With 
the present values of e 2 v 2 ~ 10~ 13 and those expected in the near future (e 2 z/ 2 ~ 10~ 10 ), 
however, one is definitely on the safe side where (asymptotic) expansion methods make 
perfect sense. 
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Appendix A. Heisenberg-Euler Analysis 

In this appendix we will check our LO results (JoTJj) and ()60|) in an independent manner 
by using the Heisenberg-Euler (HE) effective Lagrangian [HOI EDI- This is the LO in the 
derivative expansion of the effective action (jSJ) but contains all orders in the intensity. 
It has the well-known proper-time representation j3] (we follow the nice review 

oo 

^=~— 2 j-e-H "l ah / ( , -l-^(a 2 -b 2 )|, (Al) 
8tt 2 J s ltanh(ebs)tan(eas) 3 v ; /' v ; 

o 

and is exact for constant fields but remains approximately valid for photon frequencies 
small compared to the electron mass as in (fT2^ . The quantities ±ia and ±b are the 
eigenvalues of the constant matrix 

G^v = + fp,v , (A. 2) 

consisting of both background and probe field. They are related to the standard scalar 
and pseudoscalar invariants defined in and via 

a 2 - b 2 = - ^G^GT = 2^ , (A.3) 

ab = - ~G M „G^ ee & ■ (A.4) 

Note that in this appendix the invariants 5? and & denote the contribution of both 

background and fluctuation, cf. (|A.2J) . and thus are nonvanishing. The representation 

(jA.lj) contains all orders in the field G^. Low intensities allow for a weak- field expansion 

the first two orders of which are 

a 2 / 8 ^ 2 14 a 3 /256vr s 416vr ^ _,\ 

b££ = — — S? 2 + — @> 2 + — + . A.5 

m 4 ^45 45 J m s e \ 315 315 J 

It is worth to point out that exactly the same numerical coefficients appear as in (|7T| . 
For what follows it is useful to write the LO of (jA.5j) as 

5J? = \^y 2 + l - 1+ & 2 + 0(a 3 ) , (A.6) 



with the couplings r y± given by 



7+S 7{, 7 . S 4{, (A.7) 
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To approximate the polarisation tensor IP" we start from (|A.5|) . decompose into 
background F and probe / according to (jA.2|) and expand the HE action, 5S = 
J d^xSJtf, to second order in the probe field /, 

8S[FJ} = l(f aP ,^U„) . (A.8) 

The tensor Q a/3tJ-l/ is proportional to the second derivative of <5Jzf, 

^ = 4 f/fJ ■ (A. 9) 

So far we have not exploited the fact that our background consists of crossed fields. 
Note that the derivative in (jA.9j) is evaluated right at the background. It is easy to see 
that enormous simplifications arise in the crossed-field case due to the vanishing of the 
background invariants. The generic term in the HE Lagrangian is of the form y n ^ 2m 
(as odd powers of are forbidden by CP invariance) with n and m integers. Taking 
the two derivatives in (jA.9|) at / = one will always end up with (vanishing) powers of 
y and unless n and m are sufficiently small. The only surviving cases turn out to 
be the Maxwell term, {n = 1, m — 0), and the LO (jA.6|) with n = 2, m = and n = 0, 
m — 1. Thus, for crossed fields, the effective Lagrangian (|A.1J) gets truncated after the 
LO a 1 Im\. In other words, the LO describes the exact dependence on intensity for 
crossed-field background}]! In view of these considerations the tensor Q} il ' a ^ simplifies to 

Qf*>aP = 1 _F a ? F^ + 7+F Q/3 F*™ . (A. 10) 

Note that it is symmetric upon exchanging (p.v) <-» (a (3) and antisymmetric both in \i 
and v as well as a and j3. 

In momentum space the polarisation tensor Tl^ u is then given by contracting (|A.10|) 
twice with the wave vector k associated with the probe field f^, 

nr(k) = w auf3 k a k p = 7 _ w + 7+ m v . (a.h) 

To LO in k 2 this coincides with ({21} as it should, the nonvanishing eigenvalues being 
given by n ± = 7±6 2 (/c). They imply the dispersion relations, 

k 2 -l±b 2 {k) = (A.12) 

and effective metrics 

h^ = g^ + 1± T^. (A. 13) 

Introducing the index n of refraction according to (@J) the dispersion relations (jA.12|) 
and (|A.13j) become four quadratic equations for n. Demanding n = 1 for j± — singles 
out two of them. These are conveniently written in terms of the abbreviations (|29H3H) . 

n ± = i_ 7±( ^_M fc) { \/ 1 + ~ ^ " ( k ' 5 ) 2 " HH ^ " 7±k ■ 5} (A.14) 

and have the small-7± expansion, 

n± = 1 + 7± (n - k ■ s - n k /2) + L^[(n-k-s- n k ) 2 - n k /2] +... . (A.15) 

jj An analogous observation has been made for photon splitting in a plane- wave field |14j . 
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For a head-on collision of probe and background (jA.14j) yields the simple expression 

n± = - - = \ + - i + A ± , A.16 

1 - 7 ± I 1 - 7± I 

which is exact to all orders in the intensity /. Noting that 

^ = — SS' (A - 17) 
the result (jA.16|) coincides with (joTJj) and (foTIj) . 
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